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ABSTRACT 

We report a correlation between velocity offset (f3 — v/c) of strong Mg II absorption systems 
and the bolometric luminosity (Xbol) of quasars in SDSS-DR7. We find that, j3 shows a power 
law increase with -Lboi, with a slope ~ 1/4. We find that such a relation of j3 with Lboi is 
expected for outflows driven by scattering of black hole radiation by dust grains, and which 
are launched from the innermost dust survival radius. Our results indicate that a significant 
fraction of the strong Mg II absorbers, in the range of j3 = 0-0.4 may be associated with the 
quasars themselves. 

Key words: galaxies : active - quasars : general - quasars : absorption lines 



1 INTRODUCTION 

The study of Mg II absorption line systems in the spectra of quasars 
(QSOs) has been useful in detecting distant n ormal field galaxies 
situated close to the line s of sight of QSOs dBergeron & Boissel 
ll99ll;ISteideletalJll994h . Conventionally, all such absorbers with 
velocity < 5000 km s _1 relative to the background QSO are be- 
lieved to be associated to the QSO itself ('associated systems') 
while those at larger velocity offset are believed to be entirely in- 
dependent of background QSO. This general belief was questioned 
recently by the puzzling results on the abundance of strong Mg II 
absorber having equ ivalent width (W r ) more than 1.0 A : (i) by 
IProchter et alj 1 2006) where they found 2-4 time excess of strong 
Mg II absorber towards the 7 ray burst (GRB) sources relative 
to QSO sight lines (see also ISud ilovskv et al. 2007; Verga nTet alj 
l2009t iTeios et alj|2009T) . and (ii) bv liiergeron et al] feOllh . where 
they found similar excess by a factor of about 2 (3 a confidence) 
towards 45 blazar sight lines. 

These counter-intuitive results, have inspired many alterna- 
tive explanations, such as dust extinction towards QSO sight lines 
which can lower the apparent incidence rate of absorbers, or 
gravitational lensing which can increase it toward GRBs/blazars, 
but none have been found to explain th e above discrepancie s 
jPorciani et all 120071 : iMenard et all 120081 : lLawther et all l2012h . 
However the blazars, as a class, are believed to have re lativistic jet 
pointed close to our line of sight. Berg eron et al.ld201 it) speculated 
that such powerful jets in the blazars are capable of sweeping suffi- 
ciently large column densities of gas (up to 10 18 -10 20 cm~ 2 ) and 
accelerating such clouds to velocities of order ~ 0.1c, thereby pos- 
sibly accounting for the excess of Mg II absorption systems towards 
blazars in comparison with QSOs. However, such an excess in 
number of Mg II absorbers per unit redshift (dN/dz) was not con- 
firmed in the analysis of flat-sp ectrum radio quasars (FSRQs) by 
IChand & Gopal-Krishnd 12013) . though FSRQs also possess pow- 



erful jets similar to blazars, which they reconciled with the above 
hypothesis of relativistically ejected absorbing clouds, by suggest- 
ing that perhaps due to larger angle from the line of sight (unlike 
blazars with smaller angle), these accelerated clouds might not in- 
tersect the line-of-sight in the case of FSRQ s. Using a larger sam - 
ple size of 95 GRB (in cluding 12 GRB from lProchter et alj|2006h . 
ICucchiara et alj j2Q12f) did not confirm the origi nal enhancement 
found in the case of GRB bv lProchter etai] d2006l) . though a signa- 
ture of mild excess of about 1.5 times was noticed for strong Mg II 
absorption systems, albeit with only a low confidence level of 90%. 

The firm conclusion for jet based above excess still await the 
realistic numerical modelling of jet-ambient gas interaction espe- 
cially for the excess se en towards blaza rs (about a factor of 2) 
and CDQs (about 10%) jjoshi et alj|2013l) vis-a-vis normal QSOs. 
However an alternative scenario, which could b e more plausible, 
is the dust or radiation driven outflows (e.g. IScoville & Normanl 
[1991 . For instance, if there is some contribution to dN/dz of 
strong Mg II absorber from these outflows, then one will expect that 
AGN luminosity should have statistical correlation with the veloc- 
ity offset of the strong Mg II absorber relative to the background 
AGN, which is usually defined by, 



P = 



(1 + z qso ) 2 - (1 + z abs ) 



(l + 2q 



(1 



(1) 



where ft — v/c, z qso is the emission redshift of the QSO and z a b s 
is the absorption redshift of the Mg II system. 

In this letter we report a correlation between the f3 of strong 
Mg II absorbers and the bolometric luminosi ty (Lboi) of QSOs, us - 
ing the strong Mg II absorber catalogue by lLawther et alj d2012l) . 
We also propose an explanation for this correlation which draws 
upon radiation driven outflow models. In §2 we describe the sam- 
ple of strong Mg II absorbers and our selection criteria. In §3 we 
present our results and a theoretical model of radiation driven out- 
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flows. In §4 we study the fractional number counts of absorbers, 
and discuss our results in §5. 



2 DESCRIPTION OF THE SAMPLE 

We consider a sample of 10367 strong Mg II absorbers with equiva- 
lent width W r (2796 ) > 1 A belonging to 9144 QSOs, from the re- 
cent compilation by Lawther et al. (20 121) based on 105783 QSOs 
of SDSS DR7 lAbazaiian et alj|2009l ; lschneider et alj2010h . How- 
ever, the range of /3 varies with z qso , and the observed wavelength 
range of the spectrum. Therefore, in order to make the sample unbi- 
ased, firstly, we have considered a SDSS spectral range from 4000- 
9000 A which is a little narrower (by about 100 A) than the actual 
one. We then applied the following four selection filters. 

(i) We removed 773 broad absorption line (BAL) QSOs from 
our above sample to avoid any contamination in our analysis by 
BAL features which has resulted in the removal of corresponding 
931 strong Mg II absorbers. 

(ii) For all the quasars having z qso > 2.21024, the Mg II emis- 
sion line will fall above 9000 A, which is our conservative upper 
limit on wavelength of SDSS spectrum. As a result, SDSS spec- 
tra for such sources will not allow any detection of strong Mg II 
doublet falling in the redshift range between 2.21024 up to z qso . 
Therefore to avoid this observational bias, we excluded all sources 
having z qso 2.21024 from our sample, which resulted in the 
removal of 43 QSOs having 52 strong Mg II absorbers. 

(iii) Another filter was applied to avoid the observational bias 
which might result from the lower wavelength limit, viz 4000 A, in 
the SDSS spectra. In our analysis we aim to see any correlation of 
luminosity with velocity offset up to about 0.4c. However for 4000 
A considered as the conservative starting wavelength of our spectra, 
Zqso = 1.185 will be the minimum redshift, which allows us to 
detect Mg II absorber (if any) at least up to a velocity offset of 0.4c. 
Therefore, we have removed 1461 sightlines with z qso < 1.185 
having 1544 strong Mg II absorbers in their spectra. 

(iv) After applying the above mentioned redshift cuts, we are 
left with the systems with 2.21024 > z qso ^ 1.185. In these inter- 
mediate redshift systems, the /3 value can be larger than 0.4, which 
in principle may give rise to a bias of higher /3 with increasing z qso . 
Hence we also remove all the absorbers with f3 > 0.4 from the 
remaining sample which amounts to exclusion of 1523 absorbers 
along 1439 sightlines. One should note that /? = 0.4 is chosen be- 
cause if we keep /3 value less or greater than this, then the sample 
is significantly reduced. Another motivation as will be clear in the 
coming sections, is that f3 ~ 0.4 is an upper limit for the radiation 
(dust) driven outflows. 

Finally, we are left with 6317 strong Mgll systems along 5682 
QSOs in the selected redshift range. Bolometric luminos ities for the 
QSOs in SDSS DR7 are calculated in a recent study bv lShen et alj 
( 120 111) . We cross matched the QSOs in our sample, with the cata- 
logue described in Shen's paper to obtain the bolometric luminos- 
ity. We then removed two more absorbers whose QSO luminosi- 
ties were < 10 45 erg s _1 . Our final bias free sample consists of 
6315 strong Mg II systems with luminosity range 10 4 °' 5 < Lboi ^ 
10 47 8 erg s~\ with redshift range 1.185 < z qso < 2.21024, and 
with the velocity offset range of < /3c ^ 0.4c. In Figure 1, 
the blue dashed line represents the dist ribution of s t rong Mg II ab- 
sorbers in SDSS-DR7, compiled bv lLawther et alj d2012l) . and the 
black solid line is the final sample selected for this study. 
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Figure 1. Histograms with z qso of the samples of strong Mg II absorbers in 
S DSS-DR7. The blue dashed line is for 10367 strong absorbers compiled 
by Lawther e t alj j2012l) . The black solid line represents the sample used in 
this work. 

3 CORRELATION BETWEEN (3 AND L hol : SIGNATURE 
OF RADIATION DRIVEN OUTFLOW 

In order to test the dependence of f3 on luminosity, we divide the 
sample in bins of bolometric luminosity. Most of the absorbers 
(5651 out of 6315) belong to QSO sightlines having a luminosity 
range 10 46 -10 47 erg s _1 . We divide these 5651 systems into four 
bins of bolometric luminosity. We also have two more bins, one for 
Lboi < 10 46 erg s _1 , and another with Lboi > 10 47 erg s _1 , the 
first having 27 systems and the second with 637 systems. 

Consider the case of the absorbers being distributed uniformly 
in the allowed range of z a b s (which in turn is determined from the 
allowed range of /3), then the median value of f3 should be inde- 
pendent of z qso (see Appendix A for a proof). Hence, irrespective 
of the distribution of z qao in a luminosity bin, the median value of 
beta should be same in all luminosity bins. To test this hypothe- 
sis, we estimate the median, the lower 25 percentile and the upper 
25 percentile of data in each of the above mentioned six luminos- 
ity bins. We plot the median with circles, and the upper and lower 
percentiles as the end points of vertical dotted bars in Figure 2. 

Interestingly, we find that the median is not constant. The data 
shows a correlation of ft with the Lboi- The 5651 absorbers sys- 
tems with 10 46 ^ Lboi ^ 10 47 erg s _1 , which form the mainstay 
of the sample show a power law increase of f3 with Lboi, with a 
slope of ~ 1 /4. Increase of median value of j3 with the bolomet- 
ric luminosity, proves that the distribution in each bin is not uni- 
form rando m. This fact is also hinted in the evolution of dN/dz 
with z a ba dZhu & Menarj2o"l2h . Using the evolution in dN/dz for 
our sample, we evaluated the expected relation between the median 
value of P and z qso from equation JAU . We then converted it to the 
corresponding relation between /3 and luminosity by using the best 
fit relation between z qao and luminosity, a characteristic of mag- 
nitude limited survey such as SDSS. We have shown this relation 
using a dashed curve in Figure 2. Although the dashed curve does 
show some evolution, but it is clear that it cannot fully explain the 
observed /3-Lboi correlation. 

Which physical processes can give rise to non-uniformity of 
absorber distribution? The evolution in dN/ dz has been attribute d 
to the evolution in global star formation rate (Zhu&M enardl2012l) . 
although without any concrete evidence. Also, observations of in- 
tervening galaxies show a small covering fraction (^ 0.3) for 
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Figure 2. Correlation between the bolometric luminosity of the QSO and f3 
of Mg II absorber. The circles represent the median of data in a particular 
luminosity bin. The upper and lower extreme of the dotted vertical lines 
give the location of upper and lower 25 percentile, respectively. Sizes of 
the luminosity bins are shown by the horizontal bars. The solid, dotted and 
dash-dotted line represent the theoretical model discussed in §3.1. 



strong Mg II absorbers (W ^ 1 A) ^Nielsen et all2012l ; lchen et alj 
2010). Here we explore an alternate based on outflows associated 
with QSOs, which can give rise to non-uniformity of incidence of 
absorbers. As we explore in next section, the relation fit tx L^f , is 
a natural consequence of QSO radiation driven outflows. 



3.1 Absorbers as radiation driven outflows 

Radiation driven outflows have been invoked repeatedly in litera- 
ture to explain the co-evoluti on of black hole a nd bulge, to explain 
the accretion disc winds (e.g. Proga et al]|2000b and galactic winds 
(e.g. lMurrav et al.ll2005l ; ISharma et al.l201ll) . We consider here the 
radiation driven outflows, where the photons scatter the dust grains 
and impart their momentum to dust. The dust in turn is collision- 
ally coupled to the gas, and the momentum is uniformly distributed 
to the dust and gas mixture. In this scenario, the motion of dust 
and gas mixture surrounding the QSO is governed by the following 
equation, 



dv 



GM. 



dr 



(2) 



dr 4-7rr 2 c r 2 

where M. is the mass of the black bole and <3> is the dark matter 
halo potential. L uv is the integrated UV luminosity and for QSOs 
where the main emission is in high frequency bands, luminosity 
over UV and EUV ban ds is roughly half of th e bolometric lumi- 
nosity (Luv ~ 0.5Lhni)( lRis aliti & Elvis 20(3). k is the frequency 
averaged opacity for the scattering and absorption of UV photons 
by dust grains. For wavelength of photon < 0.3 /im, the k for a 
dust and gas mixtur e ranges from 200 to as large as 1000 cm 2 g _1 
jLi&Draine]|200ih . We take a value k = 500 cm 2 g~\ which 
roughly serves as an average effective value of opacity. We can in- 
tegrate equation ((2} to obtain the following expression for velocity 



2(*(r) - Hn)) 



(3) 



where ri is the launching radius of the outflow. In the case of ra- 
diation pressure on dust grains, the opacity is generally quite high 



and hence the radiation force is many times larger than the grav- 
ity, therefore the gravitational force can be neglected. At a large 
distance the velocity attains the following terminal value 

1/2 

<'x - I ] (4) 



V 4-7TC 



rb 



The base radius (rb) for launching these outflows is an important 
factor and it should be the minimum distance at which the dust 
grains can survive. Studies on dust survival yield following relation 
between the subl imation radius of th e dust grains and the luminos- 
ity of the AGN dMor & Netzedl2012h . 

r b = Rsub ~ r b ,o ( lr , 4fi Lb °' r) ■ (5) 

\ 10 4b erg s" 1 J 

The value of rt,o is 0.5 pc for graphite grains and 1.3 pc for the 
silicate grains. Substituting equation ((5} into ((4), we obtain the fol- 
lowing expression for wind terminal speed, 



0.1c 



500 cm 2 g" 1 



1/2 



10 46 erg s" 1 



1/4 



rb,o 
0.5 pc 



-1/2 



(6) 



We note that this mechanis m has previously been discu ssed in the 
context of AGN outflows bv lScoville & Normanl (1995). These au- 
thors also arrived at similar terminal speed for a radiation driven 
outflow. 

We plot this scaling to compare with the observed correlation 
of ft and Lboi in Figure 2. The dash-dotted, solid and dotted line 
in Figure 2 correspond to rb,o = 0.2, 0.5, 1.0 pc respectively. We 
find that this simple theoretical model fits the observed correlation 
pretty well, which indicates that the absorber systems are likely to 
be radiation(dust) driven outflows. 

One is then tempted to ask as to how these outflows fit in the 
unification schemes of AGN. We find that the launching radius of 
the outflows is the dust sublimation radius, which is also the in- 
ner radius of the dusty torus. Inside the torus, the UV photons are 
quickly r eprocessed into IR. Although the IR photons ca n also drive 
outflows jDorodnitsvn et al.l201 ltlsharma & Nathll2012l) . however 
the speeds would not be large, as the IR to dust scattering cross sec- 
tion is more than an order of magnitude smaller than in UV. One 
possible way to reconcile this is the following. 

Let us suppose that the outflows do not plough through the 
main body of the torus, but consist of material lifted from the outer 
surface of the torus. In that case, as the torus material is dilute 
and highly porous at its periphery, the UV photons can in princi- 
ple travel a large distance without being attenuated and impart their 
momentum to gas and dust mixture lifted from the o uter surface o f 
the torus. More specifically, in the picture presented in lElvislfeoOOh . 
the region which we are considering should take place between the 
BAL envelope and the torus. We note that, this scenario not only 
gives rise to large velocity outflows, but it may also account for 
the small fraction (< 0.1) of the QSOs which show these absorbers 
owing to the fact that the region allowed for the outflows (periphery 
of the torus) occupies a very small fraction of the viewing angle. 



4 FRACTIONAL NUMBER OF ABSORBERS 

Next we study the dependence of fractional number counts of ab- 
sorbers as a function of QSO luminosity. We define the fractional 
number count as below, 



Frac. number count = 



Number of absorbers found 
Number of QSOs searched in a bin 
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Figure 3. The fractional number count as a function of bolometric lumi- 
nosity is shown using filled diamonds, and the size of the luminosity bin 
is shown by the horizontal bar. Thin horizontal line represents the average 
value of fractional number count over the entire sample. 
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Figure 4. Ratio of absorbers with < 0.0167 to the total number of 
quasars in a particular luminosity bin is plotted against the bolometric lu- 
minosity. The horizontal bar represents the size of the luminosity bin. 



Again, we limit our analysis to the spectral region with /3 < 0.4. 
From our sample, as described in §2, we can easily estimate the 
"Number of absorbers found" in a given luminosity bin, having 
/3 < 0.4. However to find the corresponding "Number of QSOs 
searched in a bin", we also need to count those QSOs in the parent 
sample of SDSS-DR7 from which the QSOs wi th Mg II ab s orbers 
are selected. We use the parent catalogue from IShen et al.l ( 120 1 lh 
of which the sample used in this work is a subset. Therefore, we 
estimate the "Number of QSOs searched in a bin" by using non- 
BAL QSOs from Shen et al (2012) catalogue, which satisfy the 
redshift criteria 1.185 ^ z qso < 2.21, to ensure the absence of any 
observational biases (see §2). 

We plot the fractional number count as a function of luminos- 
ity (Lboi) in Figure 3. The values are shown by filled diamonds 
whose x-coordinate is the centre of each luminosity bin. We also 
show the overall average of the sample using a horizontal line, 
whose value is ~ 0.1. We find that the fraction increases steeply 
with increasing QSO luminosity and reaches a maximum roughly 
forLboi ~ 10 47 erg s _1 . For iboi > 10 47 ' 5 erg s _1 there is a mild 
decrease with luminosity, however this decrease is uncertain as in 
this bin, we have many apparently faint high redshift quasars, for 
which the signal to noise criterion removes large chunks of spectra 
and the corresponding absorbers (D. Lawther pvt. comm.). We note 
that the fractional number of absorbers has a contribution from out- 
flowing and intervening systems, which we can not separate here. 



5 DISCUSSION 

We would like to emphasize an important point in connection with 
our result. There is a general consensus in the literature, which goes 
along the line that the absorbers which have j3 < 0.0167 (v < 5000 
km s _1 ), are associated with the QSO and with j3 higher than this 
represent the intervening media. We stress here that this criterion is 
not adequate to denote the associated systems, and the true associ- 
ated systems can also have /3 > 0.0167, e.g. the QSO driven high 
velocity outflows considered here. 

To illustrate this, we plot in Figure 4, the ratio of absorbers 
with /3 < 0.0167 to the total number of absorbers in a particu- 
lar luminosity bin as a function of bolometric luminosity. One can 
clearly see that lower /3 are possible for only lower luminosity, and 



vice versa. Firstly, the figure once again confirms that the velocity 
offset /3 is correlated with luminosity, because low j3 absorbers ap- 
pear along the sightlines of low luminosity QSOs. Secondly, this 
plot, in conjunction with the correlation of j3 with Lboi, shows that 
the systems which are really 'associated ' with the QSOs are spread 
all the way from /3 = 0.0 to 0.4. 

Our results call for a study to separate out the truly associated 
(outflowing) systems and the intervening ones. Of course, one te- 
dious way to do this is to locate the intervening galaxies in each 
quasar sightline, however yet another way can be through the de- 
tailed study of line shapes and features arising from outflows and 
intervening material. We look forward to such a study in the future. 

There is another implication of the observed dependence of 
fractional number count of absorbers on QSO luminosity. If one 
considers a sample of a particular type of QSOs that covers a re- 
stricted luminosity range, then the relative number of absorbers 
may differ for different samples, and be different from the overall 
average. If we consider the right side of Figure 3, corresponding to 
Lboi > 10 46 ' 5 erg/s, there the fractional number counts are roughly 
double of the overall averag e value of 0.1. We no te that recent ob- 
servations of ~ 45 blazars (Bergeron et al.ll201 H) report an excess 
of Mg II absorbers relative to that in QSOs. We speculate here that 
this excess may also arise from the fact that the blazar sample is 
small, and it may be possible that it is biased towards higher lu- 
minosity, where the fractional number count is larger. It is possible 
that if the analysis is repeated with a larger sample of blazars then 
the excess may fade away. In fact, a similar conclusion has been 
reached for a sample of FSRQs and 7156 lobe and core dominated 
QSOs where in both cases one finds only a mild excess (loshi et al. 
120131) . In this regard we bring a recent paper by ICucchiara et alj 
d2012l) to the attention of the reader, regarding the excess seen to- 
wards GRBs, where with a large sample of GRBs the puzzle of 
Mg II incidence rate indeed disappears, and one does not find any 
excess. 

In summary, we have found a correlation between the velocity 
offset of strong Mg II absorbers and the luminosity of QSOs. The 
velocity offset (/3c) has been found to increase with the luminosity 
with a power law index ~ 1/4. We have found that radiation driven 
outflows from QSOs can give rise to such a dependence of f3 on 
Lboi- These findings lead us to conclude that a significant fraction 
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of strong Mg II absorbers (even with v > 5000 km s 1 ) along 
QSO sightlines may be the AGN driven outflows. 

We are grateful to D. Lawther for supplying the redshift path 
data. We thank an anonymous referee for insightful comments. 
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APPENDIX A: THE MEDIAN VALUE OF /? 

From equation (1), we know that /3 mainly depends on the difference of 
(1 + 2qso) and (1 + z a bs)- Here we address the question whether or not 
the median value of f} which we have computed over its fixed range varying 
from 0-0.4, depends on the value of z qso . For maximum value of f} = 0.4, 
the lower limit of z a b s , using equation (1) is given by = (C — 1) + 
Czqso , where C is a constant whose value is 0.65 for j3 = 0.4. The median 
value of absorber redshift (z™^) is the solution of following equation, 



where dN/ dz is the number of Mg II absorbers per unit redshift. If the ab- 
sorbers are distributed uniformly and the quantity dN/ dz is constant, then 
from equation (A[) we get, z™f = ((C - 1) + (C + l)z qso )/2 . Us- 
ing equation (1) we find the corresponding median value of /3, which is 
Aned ^ 0.19 , independent of z qB0 . However, the observed /3 mc( j evolves 
with luminosity (and z qso ), which indicates that absorbers are not dis- 
tributed uniformly. 




(Al) 
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